Bounds on universal quantum computation with perturbed 2d cluster states 
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Motivated by the possibility of universal quantum computation under noise perturbations, we 
compute the phase diagram of the 2d cluster state Hamiltonian in the presence of Ising terms and 
magnetic fields. Unlike in previous analysis of perturbed 2d cluster states, we find strong evidence 
of a very well defined cluster phase, separated from a polarized phase by a line of 1st and 2nd order 
transitions compatible with the 3d Ising universality class and a tricritical end point. The phase 
boundary sets an upper bound for the amount of perturbation in the system so that its ground state 
is still useful for measurement-based quantum computation purposes. Moreover, we also compute 
the local fidelity with the unperturbed 2d cluster state, showing that the actual "useful" region for 
measurement-based quantum computation extends through almost the entire cluster phase. Besides 
a classical approximation, we determine the phase diagram by combining series expansions and 
variational iPEPS methods. Our work constitutes the first analysis of the non-trivial effect of 
few-body perturbations in the 2d cluster state, which is of relevance for experimental proposals. 
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Introduction.- What defines universal computation- 
ability of a quantum computer? Or in other words, 
when can a device be used to perform universal quantum 
compuation under practical assumptions? Certainly, this 
question is of paramount importance if we aim to develop 
large-scale quantum algorithms eventually. A way to ad- 
dress this problem is through the exciting approach of 
measurement-based quantum computation (MBQC) in- 
troduced by Raussendorf and Briegel [ ] . In this setting, 
a quantum computation is performed by implementing 
local measurements on a highly-entangled quantum state, 
known as cluster state [2]. As such, MBQC is a universal 
model of quantum computation. Therefore, it is relevant 
to know how much can one perturb the cluster state be- 
fore it can no longer be used as a resource for quantum 
computation. If we see the cluster state as the ground 
state of a many-body system, then the relevant question 
is whether the ground state of the perturbed system is 
still qualified for quantum computation or not. 

Unlike in Id, the ground state of such a "cluster" 
Hamiltonian in 2d (the 2d cluster state) is a universal 
resource for MBQC. The robustness of this state can be 
naturally analyzed at zero temperature by considering 
the effect of different Hamiltonian perturbations [3-6]. 
The case of a magnetic field along the x direction was 
first considered in Ref. [3], where the model was shown 
to be self-dual, and where the whole energy spectrum was 
itself also dual (up to degeneracies) to that of the quan- 
tum compass model [7, 8], to the Xu-Moorc model [9, 10], 
and to the toric code in a transverse field [1 1 ] . The phase 
transition in all these models takes place at the self-dual 
point and it is found to be first order [11, 12]. The more 
intricate case of combined magnetic fields along x and z 
directions was later considered in Ref. [ ]. Here a line of 
lst-order transitions is detected in the x — z plane ending 



at a quantum critical point. No sign of a phase transition 
along the z-axis or close to it is observed. 

Nevertheless, it still remains a mystery the way in 
which more general perturbations affect the quantum 
computational power of the 2d cluster phase. For in- 
stance, it is unknown whether the presence of few-body 
interactions modify significantly the phase diagram of the 
2d system. This is also of practical importance, since ex- 
perimental realizations of cluster states with quantum 
simulators may also be affected by the spurious presence 
of such few-body terms [13]. 

Motivated by the above, we analyze in this letter the 
zero-temperature properties of the 2d cluster Hamilto- 
nian on a square lattice under the effect of homogeneous 
Ising perturbations in the z direction ( J z ), and also under 
an homogeneous magnetic field in the x direction (h x ). 
Unlike in the case of having magnetic fields only, we will 
see that introducing J z stabilizes a well-defined bound- 
ary for the cluster phase, which sets an upper bound for 
the amount of perturbation in the system so that the 
ground state is still useful for MBQC. More precisely, 
starting at h x — we find a line of 2nd order transi- 
tions in the 3d Ising universality class, which switches to 
a line of 1st order transitions at intermediate values of h x 
(tricritical end point), and continues until the self-dual 
point (J z = 0, h x = J) where J denotes the energy scale 
of the cluster Hamiltonian (see below). Our conclusions 
are further reinforced by calculations of the critical ex- 
ponent zv. Moreover, we also compute the local fidelity 
[15-18] with the unperturbed 2d cluster state, showing 
that the "useful" region for MBQC extends through al- 
most the entire cluster phase. Regarding the methods, 
our analysis proceeds in two steps. First, an approximate 
dual Hamiltonian describing the clusteron dynamics is 
derived which is exact for both limiting cases that only 
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one of the perturbations is present. Most importantly, for 
h x = 0, the model reduces to the well-known transverse 
field Ising model establishing the existence of a 2nd-order 
phase transition. A classical approximation on this dual 
Hamiltonian is then done to gain intuition about the per- 
turbed cluster Hamiltonian. Second, the phase diagram 
of the quantum model is computed by combining series 
expansions [23, 24] and variational iPEPS methods [20], 
similar to the approach taken in Refs. [5, 21, 22]. Our 
work also constitutes the first analysis of the non-trivial 
effect of 2-body perturbations in the 2d cluster state. 

The model.- We consider the cluster state Hamilto- 
nian on an infinite square lattice under Ising and mag- 
netic field perturbations. The cluster Hamiltonian can 
be written as a sum of mutually commuting stabiliz- 
ers = a fi,(&j e r(fi) c 'j> w here T(fi) denotes the four 
nearest-neighbour spins of lattice site fj, and the a a are 
the usual spin-1/2 Pauli matrices with a € {x,y,z}. 
In terms of these stabilizers, the cluster Hamiltonian is 
given by i?cL = ~^Eu^/i' Adding nearest-neighbour 
Ising interactions in z-direction and a magnetic field in 
x-direction, the Hamiltonian studied in this work reads 

with J z and h x being the Ising and magnetic field 
strengths, and where the sum in the Ising term runs over 
pairs of nearest-neighbours Unless said otherwise, 

we will always consider the case J = 1. 

At h x = J z = 0, the unique ground state of the above 
Hamiltonian is the 2d cluster state on a square lattice, 
which is the +1 eigenstate of all the stabilizer oper- 
ators. For this case the model is exactly solvable, and 
excitations correspond to violations of the stabilizer con- 
straints. The latter are hardcore bosons and are called 
"clusterons" (—1 eigenvalues of some operators) [ ]. 
However, a finite value of either J z or h x destroys the 
exact solvability of the model. At very large h x the sys- 
tem is clearly in a polarized phase along the x-direction, 
whereas at very large J z the system is in a symmetry- 
broken Ising phase. The nicely entangled structure of the 
2d cluster state is therefore destroyed if these perturba- 
tions are large enough, in turn leading to a macroscopic 
rearrangement of the ground-state properties. 

Hamiltonian (1) is self-dual under the exchange of J 
and h x for any value of J z . This property is readily seen 
by the observation that the self-duality for J z = is 
proven by performing a controlled-Z operation on the 
Hamiltonian (1) [5]. Consequently, as the Ising term 
commutes with controlled- Z, the self-duality holds for 
any value J z . It is therefore important to understand 
how these different limits are connected when both per- 
turbations are present and to what extent the computa- 
tionability of the perturbed cluster state is affected. 

In order to gain a first understanding of the phase di- 
agram, we switch to an approximative description of the 



clusteron dynamics which contains the energetic prop- 
erties of the two limiting cases, where only one pertur- 
bation is present, exactly. This is done by the use of 
pseudo-spins 1/2 on an effective square lattice which 
contains the eigenvalues ±1 of operators. The clus- 
ter Hamiltonian is therefore mapped to an effective field 
term. The action of the two perturbations J z and h x is 
as follows. The Ising coupling J z flips the two eigenval- 
ues of nearest-neighbor Ku operators. In contrast, the 
magnetic field h x on a site i flips the four eigenvalues Kj 
where j is a nearest neighbor of i. In total, one obtains 

H eS = -J^2T*-J z J2<T?-h x T Jl T L T 33 T H- 

1 ijaer(i) 

(2) 

For h x — 0, Eq. (2) corresponds to a transverse field 
Ising model. In the other limit J z = two decoupled Xu- 
Moore models are found. Let us stress that Hamiltonian 
(1) is isospectral to Eq. (2) for these two cases which is 
not true for the general case because the sign of matrix 
elements not only depends on the clusteron configuration. 

From this self-duality important conclusions can be 
drawn. The perturbed cluster Hamiltonian for h x = 
is isospectral to the transverse field Ising model on the 
square lattice. This implies a 2nd-order quantum phase 
transition between the cluster phase and the symmetry- 
broken Ising phase. The latter transition is therefore fun- 
damentally different to the strong lst-order phase tran- 
sition between the cluster phase and the polarized phase 
at the self-dual point J = h x for J z = 0. 

We start by a classical approximation (CA) of Hamilto- 
nian (2) which is well suited for a CA because each phase 
of the model has a classical analogue. This is different for 
the original problem (1). Here the maximally-entangled 
cluster state has no classical counterpart. In the dual 
language, the cluster phase corresponds to a polarized 
phase where all pseudospins point in z-direction. If only 
h x is finite, one has two degenerate classical ground states 
where all pseudo-spins point either in x- or — x-direction. 
Finally, there are 196 classical ground states for the case 
when only h x is finite. Indeed, one has two decoupled 
sublattices and on each sublattice 14 classical ground 
states exist (number of spins pointing in — x-direction 
should be even for each sublattice). 

Interestingly, the CA can be done with a single-site 
unit cell. This is a consequence of the fact that the ener- 
getic effect of the effective field term cx J on pseudo-spins 
pointing in x- or — x-direction is exactly the same. One 
therefore replaces (t x ,t v , t z ) by a classical vector of unit 
length (sin 9 cos 4>, sin 9 sin</>, cos 9). The classical energy 
per site eg 1 is then given by 

ep 1 = — J cos 9 — J z (sin9 cos <j)) 2 — h x {shi9 cos 0) 4 . (3) 

The angle <p is obviously zero. Minimizing the remaining 
expression, one finds a line of 2nd-order phase transition 
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FIG. 1: Second-order phase transition for J z — /isin(a), h x = 
hcos(a), with a — 1.0 and J = 1. We detect the transition 
at /i crit « 0.36 (dashed line). 



FIG. 2: First-order phase transition for J z — hsin(a), h x = 
hcos(a), with a = 0.2 and J = 1. We detect the transition 
at ft* w 0.82 (dashed line). 



emerging from the pure Ising case h x = which switches 
to a line of lst-order phase transitions for intermediate 
values of h x / J z ending at the limiting case of two decou- 
pled Xu-Moore models for J z = (see Fig. (3)). Most 
importantly, the CA captures both limiting cases qualita- 
tively correct. Moreover, a similar behaviour of the phase 
boundary is observed for the full quantum problem (1) 
which we focus on in the following. 

To this end a number of observables for the Hamilto- 
nian (1) are determined in the thermodynamic limit. In 
particular, we will consider the ground state energy per 
site eo = limff- yoo (H)o/N and the fidelity per lattice site 
with the cluster state d = limAr^oo KV'CslV'o)! 2 ^, where 
( ■ )o denotes the expectation value in the ground state 
|V>o)) N is the number of sites in the system, and |^>cs) 
is the 2d cluster state on the infinite square lattice. We 
will also be interested in the behavior of the one-particle 
gap A = (ei — eo) - where ei is the energy per site of 
the first excitation -, as well as the critical exponent zv 
defined by A ~ (J z - jf lt ) zu at fixed h x . 

Series Expansions + iPEPS.- Our approach to com- 
pute the phase diagram of the model in Eq. (1) is a 
combined analysis of series expansions and iPEPS [20], 
which has already been successful in determining phase 
diagrams of a variety of systems [5, 21, 22]. Following the 
same approach as in Ref. [5], we have computed a high- 
order series expansion of the ground state energy per site 
eo, the one-particle gap A, and the fidelity per site d. 
The series expansion is done about all three possible lim- 
its. Due to the self-duality, the large h x expansion can 
be obtained from the expansion inside the cluster phase 
for large J. The energetic properties of the system have 
been computed using a partitioning technique provided 
by Lowdin [23] , and the fidelity has been calculated using 
a projector method introduced by Takahashi [24]. From 
the behavior of the gap, we have also estimated the crit- 
ical exponent zv. Regarding the iPEPS approach, we 
have employed the simplified update method in terms 
of Projected Entangled Pair Operators (PEPO) as de- 



scribed in detail in Ref. [5]. This scheme allows us to 
implement easily the evolution in imaginary time driven 
by the 5-body, 2-body and 1-body terms in Eq. (1). 
Also, expectation values have been extracted by using 
the directional-CTM method from Ref. [ ]. Using this 
approach, we have computed eo and d [15-18]. 

The combined approach of series expansions plus 
iPEPS is done as follows. High-order series expansions 
of the one-particle gap A allows the location of 2nd- 
order quantum phase transition points. Setting J z = 
hsin(a),h x = /icos(ce) and fixing the angle a, the crit- 
ical perturbation value h mt is defined by A(/i crlt ) = 0. 
This can be determined with good accuracy by resumma- 
tion techniques like dlogPade extrapolations. However, 
series expansions restricted to one limit are not able to 
detect lst-order phase transitions. We therefore define 
the field h* for which e PEPS (/i) < e^ E (/i) with h > h* 
holds. The order of the phase transition is then assigned 
as follows: if /i crlt < h* , then we conclude that there 
is a continuous phase transition at h crlt (see Fig. (1)). 
However, if h mt > h* , then this indicates that the series 
expansion has missed a level crossing in the ground state 
observed in the variational energy computed by iPEPS, 
and therefore we conclude that there is a lst-order phase 
transition at h* (see Fig. (2)). 

In practice, the series expansions in the cluster phase 
are done at order 7 and are used for the comparison with 
iPEPS. Additionally, order 10 (order 11) has been ob- 
tained for the ground-state energy per site e sing (one- 
particle gap A Ismg ) in the limit of large J z [26]. The 
iPEPS calculations is performed with bond dimension 
2 < D < 4 and with imaginary-time evolution with time 
steps of St = 10 -4 [20]. These approximations already 
provide sufficient accuracy for our purposes, with a typi- 
cal relative error in the energy per site eo of 10~ 3 — 10~ 4 . 

Results.- In Figs. (l)-(2) we show two representative 
examples of our calculations of A and eo for two direc- 
tions in parameter space. The case of Fig. (1) corre- 
sponds to a — 1.0. Several dlogPade's of the gap A 
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FIG. 3: Phase diagram as a function of X and Y for J = 1 
obtained by the combined pCUT plus iPEPS approach. The 
phase boundaries obtained within the CA are shown as dashed 
lines. Shaded region corresponds to fidelities with the cluster 
state being larger than 0.986. Left inset shows the critical 
exponent v along the 2nd-order line in the phase diagram. 
Right inset is a zoom of the lst-order transition points close 
to the self-dual line (vertical red line). 



are taken into account (upper panel) and one detects a 



critical point at K 



crit 



0.36 before 



„iPEPS 



(ft) < 4*(h), 



which is indicative of a continuous quantum phase tran- 
sition at ft crlt . In Fig. (2), however, we observe a lst- 
order transition. The calculation for this case is done for 
a = 0.2, and we can see that the gap A closes well after 

giPEPSc 



(ft) < eg b (ft), which happens at ft* w 0.82. 



The full phase diagram is presented in Fig. (3). A 
coordinate transformation is used to map the 3d param- 
eter space (J,J z ,h x ) onto a 2d triangle which is given by 
X = (1 - J + h x )/y/2,Y = (1 — J — h x )/V2 with the 
normalization h x + J z + J = 1. Our calculations show 
that, starting from the Ising quantum critical point at 
(J«it ~ 0.3285, h x = 0) [14], there is a critical line sep- 
arating the cluster phase from a polarized phase. Along 
this line the critical exponent zv is compatible with a 
value v 0.63 (assuming z = 1) indicating a 3d Ising 
universality class ending at a 3d Ising tricritical point 
having v = 1/2 (see left inset). Here it is important to 
stress that series expansions typically slightly overesti- 
mate critical exponents. We also see in the phase dia- 
gram that at the critical point ( J z w 0.27, h x w 0.4), the 
phase transition switches from 2nd to 1st order. The lst- 
order phase boundary seems to end right at the self-dual 
point (J z = 0, h x = 1) (see right inset). This qualita- 
tive change from 2nd-order to lst-order transition is also 
reflected in the behaviour of the gap A Ismg in the Ising 
phase. Indeed, the corresponding momentum of the gap 
changes suddenly from k = (0,0) to k = (tv, 7r) when 
increasing h x / J . 

Next, we show that this phase boundary constitutes 
an upper bound for the usability of the perturbed clus- 
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FIG. 4: Fidelity per site d shown for the cases J z — (upper 
panel) and h x — (lower panel). The inset is a zoom of the 
region close to the 2nd order transition. 



ter state as a resource state for MBQC. To this end the 
fidelity per site d is shown for the limiting cases J z = 
(upper panel) and h x = (lower panel) in Fig. (4). For 
J z = the iPEPS and SE data are in excellent agree- 
ment up to the self-dual point (dashed line) . The fidelity 
per site d is found to lie over the threshold for MBQC 
d > 0.986 [19] nearly within the whole cluster phase. 
For h x = we find the fidelity rapidly decreasing down 
right at the 2nd-order phase transition. As before, the fi- 
delity is over the threshold for MBQC over large parts of 
the cluster phase. Additionally, we also show the region 
where the local fidelity lies over the threshold d > 0.986 
in the full phase diagram Fig. (3) which extends almost 
through the entire cluster phase. 

Conclusions.- We have analyzed the quantum compu- 
tational power of the perturbed 2d cluster state. This 
is done by accurately computing the phase diagram of 
the 2d cluster state Hamiltonian on a square lattice in 
the presence of an Ising interaction and a magnetic field. 
First, we derived an approximative model of the clusteron 
dynamics which is exact for both limiting cases. The clus- 
ter Hamiltonian plus Ising interactions is isospectral to 
the transverse field Ising model establishing a 2nd-order 
phase transition between the cluster and the symmetry- 
broken Ising phase. The full phase diagram of the quan- 
tum model is computed by series expansions and iPEPS. 
We have found very strong evidence of a well defined clus- 
ter phase, separated from polarized phases by a line of 
1st- and 2nd-order transitions of 3d Ising type including 
a tricritical 3d Ising end point. This line of phase transi- 
tions sets an upper bound on the amount of perturbation 
that can be introduced in the system so that the ground 
state is still a valid resource for MBQC. Moreover, we 
also computed the local fidelity with the unperturbed 2d 
cluster state, showing that the "useful" region for MBQC 
extends almost through the entire cluster phase. The re- 
sults in this work constitute the first analysis of the non- 
trivial effect of few-body perturbations in the 2d cluster 
state, which is of relevance for experimental proposals 
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using for example ultracold atoms in optical lattices [13]. 
Finally, let us suggest that it would be interesting to do 
a similar analysis with other states universal for MBQC 
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SUPPLEMENTAL MATERIAL 



Series expansions.- Let us briefly present the main series expansions in the limit J ^> h x , J z (in units of J = 1). 
We obtained the groundstate energy per site 

e (7) - -l--h 2 --J 2 -—h i ~-J 2 h 2 -—J A --J i h - ft 6 - J 2 h 4 

~ 8 x 2 z 1536 x 32 z x 32 z 64 z x 98304 x 36864 z x 

J z A h x 2 - — J/ - J, V» - — (4) 

8192 128 294912 64 w 

and the gap up to order 7: 

^ = 2 4 J z — — /ij; — 2 J z + — J z h x — — J z h x 3</ z Y28 ^ x 128 z ^ x ~8 ^ z ^ x 

- 9 j/ + ^ j z /i 2 4 - ^ - ^ j,w - ^ - 11 j z 5 - ^ ^ 

2 1152 256 512 32 12288 

197669 2 4 _ 3143 3 3 _ 72887 4 2 _ 14391 5 _ 2625 6 324187 6 

110592 * x 1536 z x 3072 2 512 2 x 128 z + 2654208 

_ 1217425 2 5 5806163 3 4 11620739 4 3 307057 5 2 155251 6 

"1327104 Jz x ~ 5308416 z x ~ 442368 z * ~ "9216" z x ~ ~W4S~ z x 

-HIZlj 7 (5) 
256 Jz ■ (5j 

Furthermore we computed the fidelity per site in order 6 as 

,( 6 ) _ 1 _lr2_l h 2_93 4 _W_ 2 2 _ J37_ 4 _ 113 4 2961 6 271261 4 2 
8 64 256 512 36864 256 2048 589824 

6007 t 2 7 4 20171 t 6 



J * ^ - 77Trc^7g ^ . ( 6 ) 



589824 14155776 

and we also obtained the groundstate energy per site in the limit J z ^> h Xl J (in units of J z = 1) in order 10: 

(10) _ 2_ _5_, 4_^^, 6 _ 905 8 _ 1388129 io_ in , r ^ 25 3 

6 ° ~ ° * 192 x 12288 x 28311552 * 1630745395200 x iU ^+ 16 ft * J 

845 5 2023 7 305 2 ,2 14155 4 2 274619 6 2 25 3 20395 3 
h x J — h x J + —— n x J — 7-—— n x J + h x J + — h x J — ^ttt "a; J 



3 7 3 



18432 1179648 96 12288 2359296 16 9216 

1330063 , 5 t3 14155 , 2 t4 28578437 , 4 t4 845 , , 1330063 , 3 , 

H ft, 5 J 3 ft, J 4 H h x 4 J 4 ft j. J 5 H ft, J 

1179648 12288 14155776 18432 21 1179648 

274619 2/6 2023 7 13495919 9 31945903 g 2 41771704639 7 3 

+ 2359296 x ~ 1179648 21 + 163074539520 * ~ 543581798400 x ~ 203843174400 x 
1035254507081 6 4 32158658557 5 5 1035254507081 4 6 41771704639 3 7 

~ 815372697600 x ~ 15099494400 x 815372697600 x ~ 203843174400 x 

31945903 2j8 13495919 9 2 5 4 5 6 905 8 

543581798400 x 163074539520 x '" 192 12288 28311552 

1388129 10 

1630745395200 ' 1 ' 

In this limit we found a jump of the critical fc-mode which corresponds to the gap. In the following we present the 
energies of the first excited state (which implies the existence of four quasi particles in the Ising-phase) for both 
fc-modes (0, 0) T and (n, n) T in order 8: 

.isine „ 3 t2 13 , T 3,2 43 t4 119 , T , 65 , , T? 119 , , T 43 , 4 

Ajn = 8 - - J 2 + — h x J - - h 2 + J ft 2 J 3 + fta; J ft^ J + ft, 4 

(0,0) 4 2 4 768 64 384 64 768 

19993 Tfi 154661 , . 564215 , , t4 806611 , 3 t3 564215 , 4 t2 

J b H h x J ft, J H ft, J fta; J 2 

884736 442368 884736 221184 884736 

154661 5 19993 6 82873487 8 153477341 7 2091266777 2 j6 

+ 442368 x ~ 884736 x + 10192158720 ~ 1274019840 x + 2548039680 x 
4309928299 3 5 1884885289 4 4 4309928299 5 3 2091266777 „ 2 

~ 1274019840 x + 1019215872 x " 1274019840 2 + 2548039680 x 
153477341 , 7 82873487 , 8 

^ J + , nin ,K, w A - (8) 



1274019840 10192158720 
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respectively 



A? ins , = 8+i J* + \h x J+\h x 2 - ^J A -f:h x J z - 2 ^-h x 2 J 2 -f-hJ i J- ±.h x 4 + 7 J— J 6 
(vt.tt) 4 2 4 768 64 384 64 768 32768 

5975 , t5 14027 , 2 t4 57169 , , tS 14027 , 4 t2 5975 , 5 T 5 , 6 

H ^ J 5 H h x ,P H /la J H + /i x J H 

147456 32768 73728 32768 147456 32768 

17737 8 146069 7 22531727 2 6 73353331 3 5 33190037 4 4 

~ 1132462080 ~ 141557760 K ~ 283115520 11 ~ 141557760 x ~ 37748736 x 
73353331 ^ j3 _ 22531727 _ 146069 17737_ ^ _ 



141557760 283115520 141557760 1132462080 

All series in this limit are clearly symmetric under the exchange h x -n- J due to the self-duality of the model. 



